Causal Dissipative Cosmology With Variable G and A 
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Abstract 

A cosmological model with variable G and A is considered in the framework 
of Israel-Stewart-Hiscock(ISH) causal theory. Power law as well as inflationary 
solutions are obtained. The gravitational constant is found to increase with time. 
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The role of dissipative effects on the evolution of the universe has been investigated 
by many works [6,8,11,15]. Weinberg has stressed that the entropy of the early universe 
may have emerged from bulk viscosity dissipation [12]. The bulk viscosity may arise from 
the conversion of massive string modes to massless modes as well as particle creation 
effects around the grand unification (GUT) era [6]. It has been shown by Murphy and 
Barrow[8,ll] that bulk viscosity can lead to inflation. In a recent paper [10], we have shown 
that this solution also arises in the presence of variable gravitational and cosmological 
constant. In the same view Singh et al [13] obtained a similar solution with different 
energy conservation. Most of these investigation are based on Eckart theory. However, it 
is shown that this theory is non-causal and all of it equilibrium states are unstable. The 
presently available theory considering dissipative effects in the universe is the full Israel- 
Stewart-Hiscock casual thermodynamics [2,3]. In this paper we would like to investigate 
the evolution of bulk viscous stress in this casual theory where both G and A vary with 
time. Our investigation shows that the variation of the gravitational constant is also 
consistent with the full casual theory. For age values consistent with present data the 
theory allows only increasing G for the present phase. We further show that the functional 
dependence of the coefficient of bulk viscosity assumed by many workers is an additional 
ad hoc assumption. We have found power law as well as inflationary solutions in the full 
casual theory. 
The Field Equations 

In a Robertson Walker(RW) flat metric, the Einstein field equations with a perfect fluid 
energy-momentum tenser and variable G and A yield 

3^ = 87rGp + A , (1) 

3^- = -4nG(3p + p) + A (2) 
ti 

and 

,G A 



3(p + p)R=-(-p + p+—)R. (3) 
The usual energy- momentum conservation, Tf^ = 0, leads to 

P + 3|(p + p) = 0. (4) 

In the presence of bulk viscous stress (II), the pressure term becomes p — > p + IT. 
Section A 

Equation (3) and (4) yields 

p + 3H(p + p) = -3UH (5) 

and 

87rGp + A = . (6) 
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where a dot denotes a derivative with respect to time t. The density p and the pressure 
p are related by the equation of the state 

p = (7 - l)p , (7) 

where 7=constant. Following [10,13] we take 

A = 3{3H 2 . (8) 

Using eq.(6), eqs.(l) & (8) yield 

G = BH-W-®, B = const. , (9) 
p = DH 2/[1 ~ P) , D = const. (10) 

while eq.(2) gives 

H + 37(1 2 ~ ® H 2 = -AnGU, (3 + 1 . (11) 
The causal evolution equation for the bulk viscous pressure is given by [2,3] 

IT + rri = -%E - e -rU(3H + 1-^-1), (12) 

where T is the temperature, ( the bulk viscosity coefficient and r the relaxation time. 
Setting e = 0, 1 gives the truncated and the full-ISH causal theory, respectively. Now 
consider the parameterization 

C = rp, and T = T oP r , (13) 
where r = T are constants. Upon using eq.(9), eq.(ll) yields 

H + g2£L~g) ff2 = _4 7rSfr 2^/(i-« n> /3 ^ 1 . (14) 

We will consider here two cases: 
Case (1): 

(a) H = H = const, and f3 ^ 1, so that i? oc exp(H t), i.e., an inflationary solution. It 
follows from eqs.(8)-(12) that 

A = cons., G = const., p = const., II = const., (15) 

and 

C = const., and r = const. (16) 

We see that the bulk viscosity remains constant during the inflationary phase. This 
solution is obtained by Wolf for the non-causal theory and by Banerjee and Beesham for 
the casual theory. 

(b) H = const, and (5 — 1: We see that A = 3H 2 , R oc exp(y^), G = const, and p = 0. 
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From eqs.(5) and (13) it follows that II = 0, ( = and T = 0,so we get the perfect fluid 
solution. Therefore, during this de Sitter solution the bulk viscosity vanishes. Also, if 
a = ^p^-, II = and we recover the model in [9]. 
Case(2): 

A power law form of R, i.e., R = At a , A, a = consts.. Applying this in eq.(15) one obtains, 

n = -Nt- 2 '^\(3^ 1 , (17) 
where N = \-*+W-ft\f + ™^\ Equations (9) and (10) give 

p = B't- 2 ^-P\G = CiW-fl C7, B' consts., /3 ^ 1 . (18) 

From the above equation we see that II = — ap, a = const. Using eqs.(12) and (17), eq.(13) 
gives 

r = Bt, B = const. (19) 

and 

c = Cor (i + /j)/(i-0 Co = constj (3 ^ l . (20) 

It is clear that the viscosity is decreasing with time. This solution is obtained by Desikan 
[5] for Brans-Dicke theory with constant deceleration parameter under the framework of 
ISH casual theory. This is manifested in the replacement of [ct/(l + b)] by —2/3/(1 — (3). 
We see that these relations are valid irrespective of the functional dependence form of ( on 
p. We would like to mention that Banerjee and Beesham followed a lengthy way to arrive 
at this result. They, however, claimed that this solution is possible only for particular 
behavior i.e., ( oc p 1 / 2 . In fact this solution is valid for all values of [q, see eqs.(2.8) and 
(2.14)] 
Section B 

Another way of looking at eq. (3) is by assuming that 

p + 3 7 tfp = (21) 

and 

I^S^- 3 ™' < 22) 
Using eqs.(l) and (8) this can be written as 

( T^f + 37iJ)p = " 3n ^ (23) 

We now assume a power law for R, i.e., R = At a ,A,a consts.. Equation (21) and (23) 
now yield 

IT = -Np, and p = Dt" 3 " 7 , D = const., (3^1, (24) 
where N = ~ 2 + a 3 ( °l ( ^~ /3) . Using eqs.(13) and (24), eq.(12) gives 

r = Bt (25) 
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where B' 1 = 3a[ 7 (l - f (1 + r)) + (± - f )]. 

In section A, we obtained a similar relation for the relaxation time and upon using eq.(13) 
we obtain 

( = rp = Vot 1 ' 3 ^, Vo = const. (26) 
From eqs.(l), (8) and (24), the gravitational parameter is given by 

G oc t~ 2+3ai (27) 

Therefore, G increases or decreases with time according to whether 3a7 > 2 or < 2. 
However, G remains constant in the pure radiation and matter epochs. We therefore, see 
that the casual bulk effects perturb the usual FRW solution for the radiation and matter 
phases. For the present era, the gravitational parameter is increasing with time( to get a 
reasonable age parameter for the universe). 
Inflationary Solution 

This is obtained when H = H = const, i.e., R = R exp(H t), R = const In this case 
eq.(20) gives 

p = D exp(-&yH t), D = const (28) 

so that 

G = Fexp{3-fH t), F = const. (29) 

and from eq.(23) 

IT = - 7 p. (30) 

Applying this to eq.(12) we get 

r = const and ( oc exp(— 37H t) . (31) 

Note that during inflation both ( and p decreases with time exponentially. We also note 
that this solution is applicable irrespective of the value of e, i.e., for e = 1, or 0. 
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